Current induced rotational torques in the skyrmion lattice phase of chiral magnets 
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In chiral magnets without inversion symmetry, the magnetic structure can form a lattice of mag- 
netic whirl lines, a two-dimensional skyrmion lattice, stabilized by spin-orbit interactions in a small 
range of temperatures and magnetic fields. The twist of the magnetization within this phase gives 
rise to an efficient coupling of macroscopic magnetic domains to spin currents. We analyze the re- 
sulting spin-transfer effects, and, in particular, focus on the current induced rotation of the magnetic 
texture by an angle. Such a rotation can arise from macroscopic temperature gradients in the sys- 
tem as has recently been shown experimentally and theoretically. Here we investigate an alternative 
mechanism, where small distortions of the skyrmion lattice and the transfer of angular momentum 
to the underlying atomic lattice play the key role. We employ the Landau-Lifshitz-Gilbert equation 
and adapt the Thiele method to derive an effective equation of motion for the rotational degree of 
freedom. We discuss the dependence of the rotation angle on the orientation of the applied magnetic 
field and the distance to the phase transition. 
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I. INTRODUCTION 

In a magnetic metal, angular momentum can be trans- 
ferred from spins of the conduction electrons to the mag- 
netization and vice versa. In non-equilibrium condi- 
tions, this flow of angular momentum - the so-called 
spin-transfer torques - results in very interesting phe- 
nomena if the magnetization is spatially non-uniformii 
For example, a spin-polarized current is able to induce 
domain- wall motion in nano wires microwave oscilla- 
tions in magnetic multilayers?^ or vortex oscillations in 
magnetic nano pillarsi^ The ability to control magnetic 
configurations by electric currents may have interesting 
applications for non-volatile magnetic memoryii 

Most experimental studies investigate spin-transfer 
torque effects in nano structures, which will be important 
for the design of future memory devices. In such nano 
structures, rather large current densities (typically larger 
than lO^^A/m^) are needed to induce, e.g., the motion 
of domain walls but they can be applied in these sys- 
tems without substantial Joule heating. Recently, spin 
transfer torque effects at much smaller current densities 
(lO^A/m^) have been observed by Jonietz et al^ with 
neutron scattering in a bulk sample of MnSi. In this ma- 
terial a peculiar magnetic structure, a lattice of magnetic 
whirls or 'skyrmions', is stabilized in a small range of 
magnetic fields and temperatures, see Fig. [1] A rota- 
tion of this skyrmion lattice by a finite angle is observed 
experimentally if the current density exceeds a critical 
threshold value. 

MnSi is an example of a chiral magnetic metal. While 
its Bravais lattice is cubic, the atomic structure (P2i3) 
has no inversion symmetry. It is therefore 'chiral', i.e. 
the atomic crystal and its mirror image do not match. 
The chirality of the crystal implies that the magneti- 
zation likes to twist in this material due to relativistic 
effects (the Dzyaloshinskii-Moriya interaction) , typically 
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Figure 1: Magnetic phase diagram of MnSi taken from 
Ref. Q|. At small magnetic fields close to the critical tem- 
perature Tc = 29. 5K a skyrmion lattice is stabilized, histori- 
cally denoted as the A-phase. The inset shows the hexagonal 
arrangements of single skyrmions. 



by forming a spiral. Close to the critical temperature and 
in the presence of a small magnetic field, a more complex 
magnetic structure is formed: a skyrmion lattice. Such 
a skyrmion lattice is in some aspects similar to the vor- 
tex lattice in a type-II superconductor. The magnetic 
structure is organized in a hexagonal lattice perpendic- 
ular to the magnetic field and is translationally invari- 
ant in parallel direction. While in a superconductor the 
phase of the order parameter winds by 27r around the core 
of each vortex where the order parameter vanishes, the 
winding of the magnetization is more complex, see Fig.[T] 
The magnetization remains always finite but its direction 
winds once around a sphere without any singular point. 
Such a configuratioir is topological stable and is called a 
'skyrmion' after the nuclear physicist Tony Skyrme, who 
showed in a pioneering work that certain configurations 
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of pion fields have the same properties as baryonsi^iiS 

Experimentally, the hexagonal magnetic lattice in 
MnSi was detected by neutron scattering The wind- 
ing of the magnetization was identified by an extra con- 
tribution to the Hall effect Also other materials with 
the same crystal symmetry show the same skyrmion 
phase, e.g. Fei_a;Coa;Sifi2rJ^ In the latter material, the 
skyrmion structure has been directly measured using 
Lorentz force microscopy by Yu et ali^ Theoretically, 
it has been pointed out in a pioneering early work of 
Bogdanov and Yablonskiii^ that Dzyaloshinskii Moriya 
interactions in chiral magnets favor magnetic skyrmion 
textures. Based on a mean-field analysis, it was, how- 
ever, arguedii that in materials with the symmetry of 
MnSi such a phase is never thermodynamically stable but 
only metastable. In contrast, we showed theoretically in 
Ref. [ll| that the skyrmion phase becomes thermody- 
namically stable in a small temperature window when 
thermal fluctuations arc properly taken into account. In- 
terestingly, the skyrmion phase becomes stable even on 
the mean field level in films with a small perpendicular 
magnetic fieldJ^ Skyrmion-like magnetic textures in chi- 
ral magnets have, for example, also been considered in 
Refs. [liH. 

From the perspective of spin-transfer torque, skyrmion 
lattice phases as in MnSi are particular interesting as the 
peculiar twist of the magnetization results in an efficient 
coupling of currents and the magnetization. Whereas for 
traditional spintronic devices such a coupling primarily 
occurs when the magnetization winds in a domain wall or 
some nanoscopic device, it extends over the macroscopic 
bulk phase for a skyrmion lattice. Indeed, the so-called 
gyrocoupling vector, introduced by Thiele^^ many years 
ago to describe the motion of magnetic domains, becomes 
proportional to the volume as the gyrocoupling per vol- 
ume can directly be identified with the skyrmion density. 
Physically, this coupling can either be visualized as aris- 
ing from Berry phases which an electron picks up when 
moving across a skyrmion texture^ or as a Magnus force 
arising from the interplay of external and circulating in- 
ternal spin currents 1^ 

The above described gyrocoupling and further dissi- 
pative forces are expected to induce a motion of the 
skyrmion lattice above a critical current strength deter- 
mined by the pinning of the magnetic structure by disor- 
der. The resulting translational motion of the magnetic 
structure is, however, extremely difficult to observe with 
neutron scattering and has not yet been detected. In- 
stead, the most pronounced effect of a current is a ro- 
tation of the magnetic lattice by an angle as described 
above. In Ref. [8| it was shown experimentally and ex- 
plained theoretically that such a rotation arises from the 
interplay of spin-torque effects and thermal gradients in 
the sample. In this paper, we want to investigate other 
forces which can also induce rotations of skyrmion lat- 
tices even in the absence of thermal gradients. Such forces 
can arise due to small distortions of the skyrmion lat- 
tice induced by the underlying atomic lattice. Rotations 



without thermal gradients have not yet been observed 
experimentally but may become important in future ex- 
periments and/or other materials with skyrmion phases. 

In the following we will first introduce our theo- 
retical framework based on the appropriate Ginzburg 
Landau theory (Sec. |TI]) and the standard Landau- 
Lifshitz-Gilbert (LLG) equation (Sec. IIIip for the 
magnctizatioui^"— Possible modifications of the LLG 
equations due to the presence of spin-orbit coupling^ 
and its ramifications are left for future studies. We then 
use the method of Thiele^^, that means we project the 
LLG equations onto the translational mode to derive an 
effective equation of motion, from which we can infer the 
drift velocity of the skyrmion lattice (neglecting the ef- 
fects of pinning by disorder). Extending this method, we 
also derive an effective equation for the rotational degree 
of freedom. In section HVl we finally apply our theory to 
the skyrmion lattice. As the relevant distortions of the 
skyrmion lattice depend sensitively on the direction of 
the magnetic field, we derive specific predictions for the 
dependence of the rotation angle on the orientation of 
fields and currents. 



II. GINZBURG-LANDAU THEORY FOR THE 
SKYRMION LATTICE IN EQUILIBRIUM 

We begin with a brief review of the Ginzburg-Landau 
theory for the skyrmion lattice used in Ref. [lf| and 
discuss additional terms that orient and distort the 
skyrmion lattice. As we will later show, these latter 
terms are necessary to enable angular momentum trans- 
fer from the magnetization to the atomic crystal lattice. 
In the following, we will always consider chiral magnets 
with the same symmetry (P2i3) as MnSi or Fei_a;Coa;Si 
where skyrmion phases have been observedii^i^ii^ and 
are predictedii to occur generically. 

As the skyrmion lattice phase occurs only in a small 
temperature window close to the classical phase tran- 
sition, one can use a Ginzburg-Landau model to de- 
scribe the equilibrium properties. The weak spin-orbit 
coupling Aso in MnSi gives rise to a clear separation 
of energy scales that allows a classification of terms in 
the Ginzburg-Landau free energy in powers of Aso- The 
strongest energy scale is determined by ferromagnetic ex- 
change interactions that favor spin alignment, while the 
relativistic rotationally invariant Dzyaloshinskii-Moriya 
spin-orbit interaction, D ~ Asoi favors chiral spin align- 
ment on a weaker scale. Keeping only terms up to order 
AgQ, the free energy is still rotationally invariant. 



F[M] = J d^r (roM^ + J{VMf 

+ 2DM -{V X M) + UM^ - B ■ M], (1) 



where M{r) is the local magnetization, B the external 
magnetic field and tq, J, D, U are parameters {U, J > 0). 
We will choose D > that selects in the helical phase 
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a left-handed spiral with wavevector Q = \Q\ = D/J. 
As all magnetic structures develop on the length scale 
1/Q ^ 1/ D, each gradient term contributes with a power 
Aso, so that Eq. ([T]) becomes indeed of order AgQ. As we 
will discuss below, the magnetic structure is only oriented 
with respect to the atomic crystal lattice by weaker terms 
that are of higher order in Aso a-nd break the rotational 
symmetryj^i^^ 

After rescaling length f = Qr, the magnetization M = 
[f//(JQ2)]i/2M and field B = [U/{JQ^f]^/^B the free 
energy functional reduces to 



However, for a small range of B and t the mean-field 
energy difference between the latter and the skyrmion 
lattice phase is tiny. For this reason, it is important to 
consider corrections beyond the mean-field approxima- 
tion. In Ref. d^l we have calculated the correction to 
the free energy arising from Gaussian thermal fiuctua- 
tions around the mean-field solution which turn out to 
stabilize the skyrmion phase. Due to the fiuctuations the 
skyrmion lattice therefore becomes a global rather than 
local minimum. 



F = 'y d^f 



(l + t)M +{VMf 



,~,4 



+2M- (V X M) + M - JB-M , (2) 

where 7 = J^Q/U, and t = ro/{JQ'^) — 1 measures 
the distance to the critical temperature (i.e. within 
the mean-field approximation, the system is spiral spin- 
ordered for t < 0, B = and paramagnetic ioi t > 0, B = 
0). From now on, wc will omit all tildes to simplify the 
notation, but keep in mind that wc have chosen particu- 
lar units. 

To describe the skyrmion lattice at finite B, one can 
expand the magnetization in plane waves 



M(r) = M 



(3) 



where the sum extends over all elements Qj of Lji, that 
denotes the reciprocal lattice without (5 = 0, and wc 
introduced the uniform ferromagnetic component Mf. 
The Fourier composition converges rapidly when more 
and more Qj are included (as has been shown explicitly 
Ref. [3^) as the skyrmion lattice is a smooth, singularity- 
free texture. For the rotationally invariant free energy 
functional Eq. ([1]), i.e. to leading order in Asoj the recip- 
rocal lattice is a two-dimensional hexagonal lattice per- 
pendicular to the external field B, and Mf is parallel 
B. Using this ansatz in Eq. ([T]), one obtains a local min- 
imum of the Ginzburg-Landau free energy for a range of 
parameters. This minimum is characterized by an almost 
constant amplitude of the magnetization and describes in 
real space a lattice of skyrmions. The skyrmion is topo- 
logically characterized by an integer winding number (for 
B\\z) 



w j dxdy h ■ (a^fi X dyti), 

uc 



(4) 



A. Orientation of the skyrmion lattice 

Within the isotropic free energy functional of Eq. ([Ij , 
the two-dimensional lattice spanned by the Qj vectors 
is always perpendicular to B but its orientation within 
this plane, described by an angle $, is not fixed due to 
the remaining rotational symmetry around the B axis. 
This rotational symmetry is however broken by terms 
of higher order in spin-orbit coupling Aso not yet in- 
cluded in Eq. ([T}, which in turn will lead to a prefer- 
ential direction of $. Due to the sixfold symmetry of 
the undistorted skyrmion lattice, $ can only be fixed by 
terms that generate an effective potential of the form 
cos(6n$ — Lpo) with n = 1, 2, ... and lpq = const. Lowest 
order perturbation theory in terms like + My + 
or (d^M)^ + (9^M)2 + [dlMf, does not produce such 
a potential. One example of a term that can lock $ in 
lowest order perturbation theory is 



Fl^IlJ d\ {{dlMf + {dlMf + {dlMf 



(5) 



The effective potential generated by this term is ex- 
tremely small because 7^ /7 ^ Aqq [note that we use the 
rcscalcd variables of Eq. ^] but there are no terms of 
lower order in Aso which can orient the skyrmion lattice. 

For the range of parameters considered in this work, 
a positive 7^ describes the experimental observatiortii 
that one of the reciprocal lattice vectors Qj tends to be 
oriented in a (110) direction (as in MnSi). A subtle prob- 
lem is the explanation of the orientation of the skyrmion 
lattice for a magnetic field in (100) direction. For this 
special case, of Eq. ^ does not pin the angle $ to 
linear order in 7^ as F^ is symmetric under a rotation 
by 7r/2 around (100) and cos[6($ 7r/2)] = -cos[6$]. 
Therefore the orientation of vectors is determined by 
effects of higher order in Aso- 



where VL = M /\M\ is the direction of the magnetization 
and we integrate over the magnetic unit cell (UC). As 
one obtains W = —1, the magnetic texture corresponds 
to a lattice of anti-skyrmions. As discussed in Ref. [33j . 
on the mean-field level there is always a solution of the 
mean-field equations with a lower energy corresponding 
to the "conical" helix with a Q vector parallel to B. 



B. Distortion of the skyrmion lattice 

Terms of higher order in spin-orbit coupling Aso will 
also distort the skyrmion lattice so that it will devi- 
ate from the perfect hexagonal structure predicted by 
Eq. ID). In lowest order in Aso? such distortions arc, for 
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example, caused by the term 
Fd^Id [ d\ [{da^Myf + {dyNPf + [d^A'Pf] (6) 



which is also written in our rescaled units. Such a term 
is consistent with the B20 crystal structure of MnSi. The 
prefactor 713 (in rcscalcd units) is again small, 7,0/7 
Ago, but expected to be much larger than 71, of Eq. (O, 
Ih/lD ~ Ago- Note, however, that this term would lock 
the orientation $ only to order 7|j <C 7l. 

The skyrmion lattice obtained in the presence of the 
term is a two-dimensional reciprocal lattice, that is a 
distorted hexagonal lattice. Its reciprocal lattice vectors 
Qj are generically not anymore perpendicular to the ex- 
ternal magnetic field B, i.e., B ■ Qj 7^ 0, but to a slightly 
changed normal vector h. Thus, h is defined by 



n • Q, = 



(7) 



Moreover it is normalized n = 1 and has a positive 
overlap with the magnetic field J5 • n > 0. 



III. SPIN TRANSFER AND ROTATIONAL 
TORQUES 

To describe the dynamics of the orientation f2(r,i) = 
M{r,t)/\M{r,t)\ of the magnetization M(r,t) in the 
presence of spin-transfer torques due to electric cur- 
rents we use the standard Landau-Lifshitz-Gilbert (LLG) 
equation)^"— 



/3 



-aO X ( dt + -Vs'sJ]Q.. 

a 



(8) 



Here Vg is an effective spin velocity parallel to the spin 
current density, ~ Mvg ^ JcPI^i with the charge cur- 
rent density j^,, the local spin polarisation p and the elec- 
tron charge e. The last two terms describe the effect of 
magnetization relaxation, leading to a Gilbert damping 
constant a and a dissipative spin transfer torque param- 
eter /?. 

The magnetization precesses in the effective magnetic 
field iJcff ~ "17 !§• Strictly speaking Eq. ^ is only 
valid for a constant amplitude of the magnetization, 
\M\ = const, and therefore one has to define carefully 
how Hcff is obtained from a Ginzburg-Landau free en- 
ergy i^[A4'] when \M\ is varying. As shown in Ref. [s^, in 
the skyrmion phase the amplitude of the magnetization 
varies only little so that the LLG equation, which docs 
not include the dynamics of the amplitude, can be used 
as a good approximation. For our numerical implemen- 
tation, we use the approximation JTeff ~ ^JiJM 
where = (TW^) is the average equilibrium magne- 
tization; other implementations H^s will only slightly 
influence our results. 

The LLG equation, Eq. ^ , can be cast into the equiv- 
alent form 



(9) 



The skyrmion lattice of an idealized system without 
anisotropics spontaneously breaks translation and ro- 
tation invariance perpendicular to the magnetic field. 
Thiele^^ suggested in 1972 to project the equations onto 
the relevant translational modes. We will use this ap- 
proach below and extend it in a straightforward way also 
to the rotational degree of freedom. Technically, the cor- 
responding equations of motion are obtained by multi- 
plying Eq. ((9)) with Gfl, where G is the generator of 
the translation or rotation mode, and integrating over a 
two-dimensional unit cell (UC) of the skyrmion crystal. 
With the help of these equations, we derive the drift ve- 
locity Vd and the effective rotation angle 5$ induced by 
the spin-current. 



A. Translational mode 

Multiplying Eq. (|9]) with the generator of the trans- 
lational mode, which is given by G^j-g^^^fl = diCl, the 
occuring integral 

/SF 
d^r {d,n) ■ — ^ (10) 

vanishes due to translational invariance. In the station- 
ary limit, where the magnetic structure drifts with a con- 
stant velocity, fl = fl{r — Vdt), one obtains^^i^i 



with 



^ JVC 







(11) 

(12a) 
(12b) 



UC 



In Ref. [2^, the vector G (up to prefactors) has been 
identified by Thiele as a "gyrocoupling vector" as it 
translates a spin current to an effective Magnus force in 
perpendicular direction. In fact, the gyrocoupling vector 
is just proportional to the winding number W of Eq. (|4]) 
per unit cell, 



G, = AnWh, 



(13) 



and points in a direction orthogonal to the two- 
dimensional skyrmion lattice, see Eq. ([7]). It is therefore 
topologically quantized . Note that in the presence of a 
distortion term like Eq. the surface normal n is not 
necessarily parallel to the applied magnetic field B. The 
effective Magnus force caused by G must be equal to a 
counter force to the electrons. Indeed, electrons which 
follow adiabatically the magnetic texture, pick up a ge- 
ometric Berry phase which results in a transverse force 
corresponding to an effective magnetic field of strength 
|G|/(47r^uc) where Ajjc is the area of the magnetic unit 
cell. The resulting topological Hall effect of the expected 
strength has been already observed in Ref. [T^j. 
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The dimensionless matrix T> is called the 'dissipative 
tensor— as it describes together with a and /3 the effects 
of dissipative forces on the moving skyrmion lattice. As 
the magnetization direction Q, only varies within the two- 
dimensional plane of the skyrmion lattice, the 3x3 ma- 
trix X> possesses a zero eigenvalue corresponding to the 
normal direction n. Within the plane of the skyrmion 
lattice, the matrix X> is diagonal in lowest order in spin- 
orbit coupling Aso due to the 6-fold symmetry of the 
skyrmion lattice. So we can approximate to lowest order 
in Aso 



D, 



VP, 



(14) 



with the projector Pij = (1 — n • n )y. To lowest order 
in the current t^s, G and 2? can be evaluated using the 
equilibrium magnetization. 

For most magnetic bulk structures G vanishes and 
therefore Vd = VsP/a. For the topological non-trivial 
skyrmion lattice phase, wc have instead a finite skyrmion 
density W = ~1 and thus a finite G. From Eq. (fTTj) we 

then obtain for the in-plane drift velocity = Pvd in 
agreement with Ref. [34| 



II /3 II 



a- /3 



v'l n X v[ 



AirW ' 



(15) 



with the in-plane spin-velocity Vg — Pvg and h is the 
normal vector given in Eq. ([7]) . The last line is obtained 
in the limit a, f3 <SC 1, which is appropriate in the limit 
of small Aso as the damping terms arise from spin-orbit 
coupling effects. The drift velocity of the skyrmion 

lattice is not parallel to the spin velocity vl due to the 
Magnus forces arising as counter forces to the topological 
Hall effect 



B. Rotational mode 

The rotational mode differs in several aspects from the 
translational one. Most importantly, weak spin-orbit in- 
teractions break rotational invariance and therefore a ro- 
tational torque due to the current can be balanced by 
a counter torque of the underlying atomic crystal lat- 
tice. Note that an infinitely large skyrmion domain for- 
mally needs an infinite time to reorient due to the fact 
that a rotation by a small angle leads to infinitely large 
time-dependent changes of the magnetization at large 
distances and to dissipation forces. In practice, domains 
are always finite and we will therefore proceed calculating 
the change of the steady-state orientation of the skyrmion 
lattice in the presence of a small current. 

In order to derive the equation of motion for the ro- 
tational degree of freedom we need the generator of ro- 
tations. Grot- The magnetic texture rotated by a finite 
angle around the axis defined by the normal vector n of 



Eq. (O is generally given by fl (r) = Rn{4))il{Rf^^ {(j))r) 
where Rfi{(f>) is the rotation matrix. For infinites- 
imal angles this rotation matrix reads {Rn{(j)))ij ~ 
6ij + (jytikjUk + 0(0^) with the Levi-Civita tensor tikj- 
From this follows the needed generator tl (r) — ^{i') ~ 
(fiGrot^ + 0{(tP) with Grot^^ = nxn~ {h(r x V))f2. As 
above, we multiply Eq. by Giot^^ and integrate over 
space to obtain an effective equation for the rotational 
degree of freedom 



Pb. (vs - Vd) + Pd il3vs - avd) = t 



(16) 



The left-hand side describes how a rotational torque is 
created by the applied current. In general, a matrix is 
needed to describe the current induced rotational motion. 
As we consider only rotations around a given axis, we can 
instead use the two vectors Pn and Pd, 



uc 



p'd = I d'rd.niG^otn), 

uc 



(17a) 
(17b) 



which wc term the reactive and dissipative rotational cou- 
pling vectors, respectively. These central quantities de- 
scribe how a velocity leads to a torque around the axis 
defined by the normal n. The reactive term arises from 
the Berry phases picked up by spin currents in the pres- 
ence of a non-trivial spin-texture while the dissipative 
terms can directly be traced back to the damping terms. 
As the magnetic texture in the skyrmion lattice phase 
only varies within a plane, the two coupling vectors are 
orthogonal to its normal, Pfl -n = and Pd - n = 0. 

As a consequence, only the in-plane velocities v^j^ and 
enter Eq. ([TBI) . 

The rotational torque exerted by the current is bal- 
anced on the right-hand side of Eq. (fTB|) by the flow of 
angular momentum from the skyrmion lattice to the un- 
derlying atomic lattice. The torque 



d^UesiGroA) 



uc 



9/ 
9$ 



^XS^ (18) 



can be expressed by the change of free energy (per mag- 
netic unit cell and divided by the magnetization) upon a 
rotation of the magnetic structure by the angle $. To ob- 
tain the correct sign note that $ describes the rotation of 
the magnetization and not of the coordinate system. In 
the linear response regime, the torque r can be expanded 
in the small deviation (5$ from the equilibrium orienta- 
tion, see discussion in section III Al The restoring force 
depends on the susceptibility, i.e., the "spring constant" 



(19) 



In Eq. (|18p . we further used that the torque van- 
ishes in equilibrium for Vs = 0. 

When discussing Eqs. ([T7|) a careful interpretation of 
the terms linear in r are necessary which arise for any 
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rotational mode. We have checked that they vanish 
for symmetrically shaped macroscopic domains. These 
terms give, however, extra rotational torques for domains 
with asymmetric shape with shape-dependent sign and 
strength. Assuming that these average to zero, we ne- 
glect all terms linear in r and approximate 

P'^ « j (fr (Ci X d.ClJ (nxCiJ , (20a) 
P"d~ [ (fr^^n(hxn). (20b) 

Within linear response, we can solve Eq. ([T5| for the 
rotational angle (5$, (provided that |x| ^ 0) 

(5$ = \Pr {v, - Vi) + Pu iPv, - avd)] , (21) 
X 

where in linear order in Vs the coefficients are again eval- 
uated with the equilibrium magnetization. Together with 
the equation for the drift velocity, Eq. (|T5|) . this is our 
central result (together with its numerical evaluation dis- 
cussed below). Note that only the in- plane velocities 
and t)|| enter Eq. (PT|) and generate a rotation as the 
orthogonal components are projected out. Using the ex- 
plicit solution for the drift velocity, Eq. (ITSI) , the equation 
for (5$ simplifies for small a and /? to 

S^^^^^f^Pj,(^nxv,)+PDVs) (22) 
X V 47r VK / 

with W = —1 for the Skymion lattice. As many other 
spin-torque effects?^ the rotation vanishes for a = (3 
where the effective Galileian invariance of Eq. ([S]) allows 
only for trivial solutions where the magnetic structure 
drifts with the current. 



IV. THEORY FOR 5$ IN SKYRMION 
LATTICES 

For the undistortcd skyrmion lattice, i.e. for the rota- 
tionally symmetric Ginzburg-Landau free energy Eq. ([l} , 
the rotation angle (5$ vanishes by symmetry. The six- 
fold rotational symmetry of the hexagonal magnetic lat- 
tice implies immediately that the two rotational coupling 
vectors Pr and P_d of Eqs. ([T7)) have to vanish. More 
precisely, as the orientation of a hexagon is described by 
a third-rank tensor, a rotational torque will only show up 
to order v^, too small for any observable effect (at least in 
bulk materials, where only relatively small current den- 
sities can be applied, sec Rcf. [H). To obtain an effect 
already in linear order inVs, one has to take into account 
that the skyrmion lattice is slightly distorted by a cou- 
pling to the underlying atomic crystal lattice. Further- 
more, one has to investigate the origin of the restoring 
forces, x^^j which also arise from higher-order spin-orbit 
coupling terms. 

We have systematically investigated which symmetry- 
allowed terms in the Ginzburg-Landau theory to leading 



order in Ago give rise to (i) a distortion allowing for ro- 
tational coupling or (ii) a preferred orientation of the 
magnetic lattice. Two representatives of such terms are 
given by Fl in Eq. (O and Fd in Eq. © with coupling 
constants 7^ and jd, repectively. F]j leads to small but 
finite rotational coupling vectors, such that Pji and P^ 
arc of order 0(715). The term F^, on the other hand, 
orients the skyrmion lattice and gives rise to a finite sus- 
ceptibility X ^ C'(7l). While 70 is very small, a sizable 
(5$ is only obtained because the susceptibility x is even 
smaller jl/jd ^ ^so- ^ consequence, a sizable effect 
(5<i> oc 70 /7l can be expected. 



A. Numerical solution of the skyrmion lattice 

For the numerical evaluation of (5$, we have to eval- 
uate the coefficients of the formula for the rotation an- 
gle, Eq. (HI]), with the equilibrium magnetization texture 
M{r) of the skyrmion lattice obtained from the free en- 
ergy functional of Eq. ^ together with Eqs. ^ and 
In order to obtain Af (r), we employ the following 
mean-field approximation. We minimize the free energy 
functional with the ansatz of Eq. ([3]) for the magnetiza- 
tion but we include in the sum only the three smallest 
reciprocal lattice vectors Qj . This is a good approxima- 
tion as it turns out that higher order terms contribute 
(both experimentally and theoreticall}>ii) maximally a 
few percent to the total magnetization. For simplicity, 
we do not include the effects of thermal fluctuations in 
our calculation. While these are essential to promote the 
skyrmion solution from a local minimum to a global min- 
imum of the free energy, they are expected to give rise 
only to small renormalization of prefactors at least not 
too close to the phase transition as was checked explicitly 
in Rcf. m. 

Hence, we will approximate the static magnetic texture 
in total by 27 real parameters fXi, i = 1, 27: the uni- 
form magnetization Af/, two reciprocal lattice vectors 
Qii Q2 (with Q3 = Qi — Q2) and three complex vec- 
tors Mq ., j ~ 1,2, 3. In this representation, the solution 
that accounts for the sliding motion of the skyrmion lat- 
tice can then be written in the absence of pinning forces 
by disorder and to linear order in the applied current as 

Mir,t)^ Mir-Vdt,{^l,}) (23) 

where is the drift velocity of the magnetic structure. 



B. Numerical evaluation of 5$ 

With the help of the magnetic texture M{r) in equi- 
librium we can evaluate the gyrocoupling vector G, the 
dissipative tensor "D, the rotational coupling vectors Pd 
and Pl, and the susceptibility x that are needed to de- 
termine the rotational angle using Eq. (|2ip . In the 
limit of smaU a, /3 < 1, Eq. ((HI) reduces to Eq. (gH), and 
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Figure 2: Effective rotation angle (5$ in units of $0, Eq. ((25}, 
for a magnetic field, B, in [1-10] direction as a function of the 
current direction perpendicular to [1-10] for various values 
of a, j3 given in the inset. The other chosen parameters are 
t = -0.8, \B\ = 0.5^^, 7D = 0.01 and 7l = 0.001. 



after collecting all prefactors arising from the rescaling of 
variables, we then obtain 

5<i> = ^o5ip{t,B,j) for a,/3,7L,7D < 1, (24) 
where $0 is given by 



)2j3/2 



7L 



(25) 



The dimcnsionless function Sip{t, B , j) depends only on 
the dimensionless distance from the critical point t, the 
direction and strength of the dimensionless magnetic field 
B, sec Eq. ([2]), and the orientation of the current. 

We will first discuss the dependence of (5$ on the rela- 
tive and absolute orientation of the magnetic field B and 
applied current v^- Afterwards we discuss the depen- 
dence of (5$ on the distance, t, to the phase transition. 



1. Orientational dependence of 

Within our theory, the rotation angle 5$ is propor- 
tional to the product of Vs and a vector in the skyrmion 
lattice plane, which is almost perpendicular to the mag- 
netic field B with deviations of order ju. As a conse- 
quence, when rotating the current around the direction 
of a magnetic field, a simple cosine dependence is ob- 
tained. This is demonstrated in Fig. [5] where a numerical 
evaluation of (5$ using Eq. (PT|) as a function of Vg for 
various values of a and l3 is shown. The blue solid line 
corresponds to the limit of small a, /3 where (5$ reduces 
to Eq. ((M)) . 

The dependence of 6^ on the direction of the field, B, 
on the other hand, is substantially more complicated. 
The three figures. Fig. |31 Fig. S] and Fig. [51 give an 



Figure 3: Effective rotation angle 5^ in units of "I>o, Eq. (|25|l . 
for a current, Vs, in the [110] direction as a function of the di- 
rection of the magnetic field, B, perpendicular to [110] (other 
parameters as in Fig. Note that for the high symmetry 
directions (111) and (100) (5$ vanishes. Close to the (100) 
direction the effect is maximal (see text). For certain direc- 
tions, the effect depends sensitively on the size of a (in the 
reseated units P is less important) as it affects the direction 
of the drifting skyrmion lattice. 



overview on how (5$ depends on the orientation of B for 
fixed current direction. Several of the main features in 
these figures can be understood from symmetry consid- 
erations as explained in the following. Special properties 
can not only be expected when the field is oriented along 
either a two-fold (100) axis or a three-fold (111) axis, but 
also for a field perpendicular to a (100) axis as in this case 
the product of time-reversal and it rotations around (100) 
maps B (and the skyrmion lattice) upon itself. 

Fig. [3] shows how the rotation angle changes as the field 
B is rotated around the [110] axis for a current parallel to 
[110]. For this geometry, one of the reciprocal lattice vec- 
tors is in the absence of a current always oriented in the 
[1-10] direction for the parameters chosen in Fig. [3] Note 
that the directional dependence is not very universal and 
depends on the structure, sign and size of the dominant 
anisotropy terms and even the size of the damping con- 
stants. A universal feature is, however, that vanishes 
(to linear order in Vg) for magnetic fields in one of the 
two relevant high-symmetry direction (note that (110) is 
not a symmetry axis of the B20 structure of MnSi). As 
(100) is a two-fold screw axis and (111) a three-fold sym- 
metry axis, both Pjj and Pji and, consequently, also 5$ 
are zero for a magnetic field oriented in these directions 
(as Pd and Pr are perpendicular to B because n be- 
comes parallel to B for these high-symmetry directions). 
Away from the two high-symmetry directions, (100) or 
(111), we obtain finite coupling vectors, Pjj and Pr, 
and a transfer of angular momentum from the magnetic 
texture to the crystal lattice resulting in a finite S^. 

As we have limited our analysis to effects that are 
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linear in the current, a reversal of the current direc- 
tion always leads to a sign change of the rotation angle, 
(5$ ^ — (5$. To linear order in the damping coefficients 
a and /3, where Eq. (|2ip reduces to Eq. (P^ . also a re- 
versal of the magnetic field direction, B —B leads 
to a sign reversal (5$ —5$ (blue solid line in Fig. [3]). 
This is not the case if contribution of higher order in 
a, /? are incorporated (time reversal is broken not only 
by B but also by the applied current and dissipative and 
reactive forces have opposite signatures under time re- 
versal). Even for realistic values of a ~ 0.1,— the latter 
can have a large effect and can even change the sign of 
(5<i> for certain crystallographic directions, see Fig. [3l As 
a consequence, there is generically no specific symmetry 
with respect to a reversal of the current and the magnetic 
field, 5^{B,j) ^ 5$(-B, -j). 

A geometry, where the field is oriented along the [1 — 10] 
direction and the current along the [110] direction, is an 
exception as (100) is a two-fold rotation axis which al- 
lows to map B B and j —j. Precisely this 
geometry has been studied experimentally in Ref. Q, 
where the same scattering pattern was observed when 
both field and current were reversed. Within our con- 
ventions for i5<i> (defined relative to the field orientation), 
this corresponds to a reversal of 5$ when both current- 
and field direction are reversed. In contrast, our sym- 
metry analysis above shows the opposite behavior (for 
all symmetry-allowed anisotropics and even beyond lin- 
ear response theory). This discrepancy was resolved in 
Ref. [8| where it was shown that additional symmetry 
breaking temperature gradients explain the experiments, 
see introduction. 

The specific directions (100) for the magnetic field 
require extra consideration as for this orientation the 
susceptibility x of Eq. (|19p vanishes to linear order in 
7l as discussed below Eq. ([5]). More precisely, denot- 
ing by 6 the angle between B and [001] (for the ge- 
ometry of Fig. [3]), the relevant potential is proportional 
to 7l(5^cos6$ while Pn and Pd vanish linearly in 5. 
Therefore (5<I> ~ 1/(5 as can be seen in Fig. [3l Only for 
very small 5 effects of order 7^ lead to a rounding of the 
divergence. liB is precisely oriented in [001], one also has 
to include further anisotropy terms in the analysis, see 
the discussion in section III Al and the predictive power 
of our theory along such special symmetry directions is 
limited. 

In Fig. |4] the current-induced rotation angle is shown 
for a current in the [111] direction when the magnetic 
field is rotated perpendicular to this direction. By sym- 
metry the pattern repeats itself every 120°. Note that 
i54> becomes small but does not vanish for fields in the 
(110) directions, see also Fig. O 

For current in the [100] direction and magnetic field 
perpendicular to [100] one again obtains large values for 
i5$ when the field points in a (001) direction. Remark- 
ably, i5$/<I>o vanishes exactly in the limit a,/3 — >■ for 
this configuration while it is finite for other orientations 
of B and j. The reason is again a special symmetry: the 
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Figure 4: Elfective rotation angle 5^ in units of "I>o, Eq. pSjl . 
for a current, Vs, in the [111] direction as a function of the di- 
rection of the magnetic field, B, perpendicular to [111] (other 
parameters as in Fig. [2]). 
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Figure 5: Effective rotation angle (5"1> in units of "^o, Eq. (|25|l . 
for current, Vs, in the [100] direction as the function of the di- 
rection of the magnetic field, B, perpendicular to [100] (other 
parameters as in Fig. [2]). Note that for this orientation, the 
rotation angle vanishes quadratically in the limit a, /3 — >■ 
due to a special symmetry, see text. 



product of time reversal T and a rotation by 180° around 
the [100] direction. Under this symmetry, B is mapped 
upon itself. It also enforces that the reactive rotational 
coupling vector, P^, which is even under T, points in 
[100] direction while the dissipative rotational coupling 
vector, Pd, which is odd under T, has to be perpendic- 
ular to [100]. Using Eqs. ([15]) one obtains that Vg — Vd 
becomes perpendicular to Vg while j3vs — avd becomes 
parallel to to leading order in f3 and a. Therefore all 
current induced torques vanish for a, /? — )■ according to 
(HSl). 
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Figure 6: Angle 5^ in the limit of small a,/3, see Eq. p4t . 
for a current in [110] direction and various orientations of the 
magnetic field, see legend, as a function of the dimensionless 
distance t from the critical point {t — —0.8, |B| — 0.5y/—2t, 
7o = 0.01 and -jl = 0.001). 



2. Dependence of 5$ on the distance to the phase 
transition and order of magnitude estimate 

When discussing the physics close to the phase tran- 
sition, one first has to emphasize that all phase transi- 
tions are expected to be of first order. This is obvious 
already on the mean field level both for the transition 
to the conical phase (where the ordering vectors jump) 
and for the transition to the paramagnetic phase (due 
to the presence of a cubic term in the Ginzburg Landau 
description in the presence of a finite magnetic field). 
The latter transition, is, however, also strongly affected 
by thermal fiuctuations which, drive even the transition 
from the paramagnetic to the helical phase first order <^ 
This effect makes it very difficult to estimate the precise 
location of the phase transition line. 

As all transitions are first order, the rotation angle is 
formally non-singular at the transition. In practice, how- 
ever, a complex interplay of phase transition dynamics, 
pinning effects, the external drive by currents, heating 
effects and even surface properties can be expected in 
the regime where both phases are locally stable. These 
questions are certainly far beyond the goal of the present 
study. 

We therefore restrict our analysis to the overall depen- 
dence of the rotation angle on the parameter t, which can 
be controlled by temperature in the experiments. Fig. [5] 
displays the dependence of (5$ on t for fixed magnetic 
field strength and various orientation of B in the limit of 
small a and /3 when Eq. is valid. For the chosen set 
of parameters, (5$ grows upon increasing t, i.e. for grow- 
ing temperature. The overall magnitude of the effect, 
might, however, be overestimated as the mean-field the- 
ory cannot describe the first-order transitions quantita- 
tively. Qualitatively different results for the t dependence 
are obtained in cases where anisotropy terms dominate 



which are not quadratic in the order parameter (as as- 
sumed by us) but are of higher order implying a larger 
sensitivity to the distance to the phase transition. 

Finally, we would like to estimate the order of magni- 
tude of the rotation. As both 7^ and jd are unknown, 
it is not possible to predict quantitatively the size of the 
expected rotation angle 5$. One can, however, make a 
crude order of magnitude estimate by counting powers 
of spin-orbit coupling Aso- For t ^ 1 we can estimate 
(5$ ^ ([)g^ see Eq. ([25|) . We approximate Q = Aso/a, 



U/J'^/'^ ^ l/ksTc, ^l/id ~ Ago (see above) and 



with the drift velocity of charge, Vs 



to obtain 



(5$ 



^so 



(26) 



For MnSi spin torque effects were observed^ for cur- 
rents of the order of 10^ A/m^. In the appropriate dimen- 
sionless units this corresponds to ~ 10~^ (using the 
lattice constant a w 4.2A and Tc ~ 30K), which shows 
the smallness of the applied currents (in most spin-torque 
experiments currents are 5-6 orders of magnitude larger) . 
In MnSi a ^ 0.1 appears to be surprisingly large as 
electron-spin-resonance experiments show a rather broad 
peak.— Using Aso ~ 0.01, one could in principle ob- 
tain sizable rotation angles ^ 0(1). Experimentally, 
no such rotation was observed in an experimental setup 
which avoids temperature gradients. Taking the crudc- 
ness of the estimates given above into account, this result 
is unfortunately also consistent with our analysis, espe- 
cially as for the experimental setup (field along [110], cur- 
rent along [1-10]) the effect turns out to be suppressed 
by another factor of 0.05, see Fig. [5] 



V. SUMMARY AND DISCUSSION 

The focus of this paper is the investigation of a specific 
mechanism how spin transfer torques can lead to a spa- 
tial rotation of the magnetic skyrmion texture by a finite 
angle. Our analysis started from the observation that 
a perfect skyrmion lattice has neither a preferred orien- 
tation perpendicular to the applied magnetic field nor 
can a small current exert (to linear order) a rotational 
torque to such a symmetric structure due to its sixfold 
rotational symmetry. The magnetic texture is, however, 
embedded in the atomic crystal of the host material. This 
environment breaks rotational symmetry leading both to 
a preferred orientation of the skyrmion lattice in equi- 
librium by tiny spin-orbit coupling effects and also to a 
small distortion of the skyrmion lattice which enables the 
current to exert rotational torques. The balance of these 
two effects determines the rotation angle. A systematic 
analysis of such effects is possible as all relevant phenom- 
ena are controlled by weak spin-orbit coupling effects and 
occur on long length scales. Overall the following picture 
emerges: First, angular momentum is transfered from the 
spins of the conduction electrons to the magnetization. 
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This induces a rotation until all rotational torques are 
balanced and the angular momentum is flowing from the 
magnetic texture (via spin-orbit coupling effects) to the 
underlying atomic lattice. 

There are several other effects which can also lead to 
rotational torques. Most relevant for the experiment in 
Ref. [sj is that macroscopic inhomogeneities can lead to 
inhomogencous forces and therefore also to rotational 
forces. In Ref. @ these could be controlled experi- 
mentally by small temperature gradients in the sample. 
A temperature gradient leads to a different strength of 
forces at the 'hot' and the 'cold' side of a magnetic do- 
main and therefore to torques. Reversing the tempera- 
ture gradient therefore leads to a reversal of the rotation 
angle. Moreover the shape of a magnetic domain can 
be the origin of forces which orient the domain in the 
presence of a current. These torques arc likely to be of 
random sign and might be responsible for the smearing of 
the neutron scattering signal observed in Ref. [1] . Finally, 
also distortions of the skyrmion lattice by disorder can 
lead to a reorientation of a sliding lattice. This physics, 
which is based on the a non-linear response of the mov- 
ing lattice, has previously been investigated for vortex 
lattices in superconductors<^2^ Finally, also the current 
itself can distort the skyrmion lattice and lead to a reori- 
entation of the lattice. By symmetry, this effect occurs, 
however, only to third order in the current density. 

Besides the more widely studied translational motion, 
we expect that also the rotation of magnetic textures 
will continue to be an important signature of spin-torque 
effects.— Depending on the setup both rotation and trans- 
lation can define 'soft modes' where pinning effects are 
weak and small forces can lead to sizable effect. To ana- 
lyze such rotational forces we have projected the widely 
used Landau Lifshitz Gilbert equation onto the rotational 
degree of freedom using a straightforward generalization 
of the approach used by Thiele^ for translational mo- 
tion. We expect that this approach should also be useful 
to analyze other sources of rotational torques both in 
skyrmion lattices and other experiments where rotation 
plays a rolci^ 

After this work was completed, a new type of damping 
term was suggested by Zang et al. in Ref. (40| which 
complements the Gilbert damping in the LLG equation. 
This additional damping will modify prefactors in our 
results for the rotation angle but we do not expect any 
qualitative changes. 
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Appendix A: Alternative derivation of the rotation 
angle 

The derivation of Eq. (pij) for the rotation angle was 
based on the projection of the LLG equation on the rota- 
tional mode around an axis defined by the normal vector 
h of the skyrmion lattice, sec [71 In order to check this 
approach and confirm the validity of Eq. (j2ip . we used 
an alternative derivation without projection to the rota- 
tional mode. 

In this alternative approach, we have solved for a 
steady-state solution of the LLG equation directly 
within the variational ansatz of Eq. (j23|) with variational 
parameters /ii, i = 1,..,A^ with N = 27. In order to 
determine the change of the variational parameters 
in the presence of a current v^, we multiply Eq. by 
dti/dfii and integrate over space. As the effective mag- 
netic field Hcff vanishes in equilibrium, we expand it to 
linear order in the deviations S/ii to obtain N equations 
generalizing Eq. ([21]) 

PrA^s - Vd) + PdAPvs - avd) = f'/jSnj (Al) 

where the generalized reactive and dissipative coupling 
vectors, P r/d,! ^-nd the stiffness matrix J^" are given by 



{PB„^)n 
{PD,i)n 

f" 

J VI 



uc 



an 

dm' 



(frdnS}, 

uc C'/^i 

uc dfii dnj 



M 



uc 



(A2a) 



(A2b) 
(A2c) 



As discussed in section [III Bl terms in the integrand that 
arc linear in the coordinate r arise from derivatives with 
respect to the reciprocal lattice vectors Qj . We again ne- 
glect such terms which is justified for symmetric bound- 
ary conditions. 

The rotation angle (5$ can be obtained from the 
changes in the magnetic texture parameterized by the 
deviations of the variational parameters. The rota- 
tion is obtained from the change SQj^ of the reciprocal 
lattice vectors of the skyrmion lattice. For small 5$ we 
have (5Qj = {dQ Jd^) 6^. Multiplying this formula by 
(dQi/d^) and summing over i allows to solve for 5$. 
Thus, we get the alternative expression for the rotation 
angle 



(5$ 



d<S> 



d<s> ' a* 



(A3) 



The deviations 6Q,i are obtained by solving Eq. (|Aip us- 
ing Eq. (|15p for the drift velocity Vd- 
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We have checked both analytically and numerically thus confirmed the validity of the Thiele approach in the 
that Eq. (|A3p and Eq. ([?!]) give identical results, and present context for the rotational motion. 
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